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Pioneering experiments on the spin structure of the nucleon performed in the sev- 
enties [1] were followed by the EMC experiment at CERN which obtained a stunning 
conclusion on the quark contribution to the proton spin [2]. This result triggered ex- 
tensive studies of the spin structure of the nucleon in lepton scattering experiments at 
CERN by the SMC [3] and COMPASS SLAC [5], DESY (HERMES) P and JLAB 
[7j as well as in polarised proton-proton collisions at RHIC [8j [9]. The results of these 
studies confirmed the validity of the Bjorken sum rule and the violation of the Ellis- Jaffe 
sum rule. In addition, the parton helicity distributions in the nucleon were extracted us- 
ing QCD analyses. The quark contribution to the proton helicity is now confirmed to 
be around 0.3, smaller than the expected value of 0.6 [10]. However, due to the limited 
range in Q 2 covered by the experiments at fixed XBj the QCD analyses (e.g. [1]) show 
limited sensitivity to the gluon helicity distribution, Ag(x), and its first moment, AG. 
The determination of Ag(x) has therefore to be complemented by direct measurements 
in dedicated experiments. 

The gluon polarisation (Ag/g) x has been determined from the photon-gluon fusion 
(PGF) process by HERMES [TT], SMC [12] and COMPASS [13]. These analyses used 
events containing hadron pairs with high transverse momenta, p<r, with respect to the 
virtual photon direction. This method provides good statistical precision but relies on 
Monte Carlo generators simulating QCD processes. The measurements point towards a 
small value of the gluon polarisation at x m 0.1. This is in line with recent results from 
PHENIX and STAR [H] at RHIC. 

In the Quark Parton Model the nucleon spin is given by the quark spins, AS, while 
AG vanishes. Taking into account orbital angular momenta, L, of quarks and gluons the 
nucleon spin is 

± = ±AZ + AG + L Z . (1) 

In QCD the U(l) anomaly generates a gluonic contribution to the measured singlet axial 
coupling, ao(Q 2 ). This anomalous gluonic contribution does not vanish at Q 2 — > oo. As a 
result, AY,(Q 2 ) becomes scheme dependent and may differ from the observable ao while 
AG is scheme-independent at least up to the NLO. In the Adler-Bardeen factorization 
scheme [H] AS AB is independent of Q 2 . Restoring the Ellis- Jaffe value of AS AB »s 0.6 
requires a value of AG(Q 2 ) m 2 and L z ps — 2 at Q 2 = 5 (GeV/c) 2 . 

Here, we present a new result on (Ag/g) x from muon-deuteron scattering. The gluon 
polarisation is determined assuming that open charm production is dominated by the PGF 
mechanism yielding a cc pair which fragments mainly into D mesons. This method has 
the advantage that in lowest order there are no other contributions to the cross section; 
however, it is statistically limited. In our analysis only one charmed meson is required 
in every event. This meson is selected through its decay in one of the two channels: 
£>*(2010) + -> D°n+i ow -> K-tt + tt+ 1ow (D* sample) and D° -> K~tt + (D° sample) and 
their charge conjugates. 

The data were collected during 2002 to 2004 with the 160 GeV/c CERN SPS //+ 
beam and correspond to an integrated luminosity of 1.7 fb™ 1 . The beam muons coming 
from the 7r + and K + decays are naturally polarised with an average polarisation, P^, of 
about -80%. 

The polarised 6 LiD target consists of two cells (upstream u, downstream d), each 

research funds for 2005-2007 
") Supported by KBN grant nr 134/E-365/SPUB-M/CERN/P-03/DZ299/2000 



1 



60 cm long, longitudinally polarised with opposite orientations. The spin directions are 
reversed every eight hours by rotating the field of the target magnet system. The average 
target polarisations, P t , were ±50%. The 6 Li nucleus basically consists of an 4 He core 
plus a deuteron. A dilution factor, /, of about 0.4 is obtained for the target material. 
The exact value is kinematics dependent and is calculated as described in |15j . Particle 
tracking and identification are performed in a two-stage spectrometer [16] . 

The present analysis selects events with an incoming muon, a scattered muon iden- 
tified behind hadron absorbers, an interaction vertex in the target and at least two addi- 
tional charged tracks. The D° mesons are reconstructed through their decay which has 
a branching ratio of 3.8%. Due to multiple Coulomb scattering of the charged particles in 
the solid state target the spatial resolution of the vertex reconstruction is not sufficient to 
separate the D° production and decay vertices. Therefore, the mesons are reconstructed 
on a combinatorial basis, considering all pairs of oppositely charged tracks in a given 
event and calculating their invariant mass. This method results in a high combinatorial 
background which is reduced in further analysis steps. 

The largest background reduction stems from kaon identification in the Ring Imag- 
ing CHerenkov counter (RICH). This restricts the studied events to a sample with kaons 
of momenta exceeding 9.1 GeV/c. Simulations have shown that in the acceptance about 
70% of kaons coming from D° decays exceed this threshold. 

Particle identification in the RICH starts from reconstructed tracks with measured 
momenta. The angles between the track and the detected Cherenkov photons are cal- 
culated for each track. The comparison of the expected angular distribution of photons 
for a pion, a kaon or a proton with the measured one is used for particle identification. 
For this comparison two different methods were used. The data from 2002 and 2003 were 
analysed using a x 2 calculated from the photon angles and the Cherenkov angles. The 
mass hypothesis with the smallest x 2 is selected. For the 2004 data, the likelihoods for the 
various mass hypotheses are also compared to the background likelihood. The background 
likelihood function is evaluated using photons not associated to reconstructed tracks. 

In the analysis an identified kaon and an identified pion are required for each event 
except for the D* sample analysed with the x 2 method where all tracks not identified 
as kaons are considered as pion candidates. The D* and the D° samples are analysed 
independently. The following two kinematic cuts are used for the D* (D°) sample: a cut 
z > 0.2 (z > 0.25), where z is the fraction of the energy of the virtual photon carried by 
the D° meson candidate, and a cut Icos^* | < 0.85 (|cos#*| < 0.5), where 9* is the decay 
angle in the D° c.m.s. system. In the D* channel a cut on the mass difference is imposed, 
3.1 MeV/c 2 < M KnWslow - M Kn — M n < 9.1 MeV/c 2 , where M KwWslow and M KlT are the 
masses of the D* and the D° candidates, respectively. The resulting signal-to-background 
ratio for events in the signal region is approximately 1:1 for the D* and 1:10 for the D° 
sample (Fig. [TJ. Note that the events entering the D* sample are not used in the D° 
sample. 

The number of events, N, collected in a given target cell and time interval (about 
one week of data taking) is 

N = a<t>n(as + <tb) x 

l + wfoLL-^M^) + D-^A B )] . (2) 
V OS + cr B \ g I x cr s + Ob / . 

Here, a, (f) and n are the spectrometer acceptance, the integrated incident muon flux and 
the number of target nucleons, respectively. In addition as (o"b) is the cross section of the 
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Figure 1: Invariant mass distributions of the Kir pairs tagged with the D* decay (D* sample, 
lower) and of the non-tagged Ktt pairs (D° sample, upper). The curves represent fits with 
functional form described in the text. 

events described by the central Gaussians (background) in Fig. 1. The analysing power 
All is the asymmetry for the fig — > ficc process, and A B the background asymmetry. The 
depolarisation factor, D, from Ref. [T7] is used. 

A simultaneous extraction of (Ag/g) x and As, assumed to be constant in the mass 
range considered, is performed independently for the D* and the D° sample using the 
events in the mass range —400 MeV/c 2 < Mpo — M Kn < 400 MeV/c 2 recorded in the 
two target cells before (u,d) and after (u',d') target spin reversal. The events of the four 
samples are weighted with a signal weight, w$, and independently with a background 
weight, w B , 

ws = Pbfa L L : , w B = P h fD ■ . (3) 

o"s + o"b cr s + a B 

The target polarisation is not included into the weights because it is time dependent. In 

this way 8 equations 




+ 7Mb (4) 



with 



Pc — ^N t l j 7c — t 
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Figure 2: Correlation between generated and reconstructed values of oll- 



are obtained from Eq. (f2]) for the 10 unknowns which are (Ag/g) x , A-q and 8 acceptance 
factors Oq = J a*0*n*(crs + a B )w t c dX with t = u,d,u',d' and C = S,B |T8] . Here, J dX 
stands for the integration over the accessible kinematic region. Assuming that possible 
acceptance variations affect the upstream and downstream cells in the same way, i.e. 
a c/ a c = a c/ a ci provides two additional equations (as in Refs. [HE]). With an extra, 
much weaker, assumption that signal and background events on the same target cell are 
affected in the same way by the acceptance variations one finally arrives at a system of 
8 equations with 7 unknowns. Possible deviations from the above assumptions generate 
false asymmetries which are included in the systematic error. 

For the evaluation of Eq. (j4]) P^ is parameterized as a function of the beam momen- 
tum. For P t , values averaged over about one hour are used. The signal purity crs/(°~s + &b), 
is obtained from a fit to the — Mpo spectra. The fit is done separately for the events 
originating from the two target cells as well as for five (three) separate bins of /Pb^LL 
for the D* (D°) sample. This takes care of the correlation between the signal purity 
and the analysing power. The spectra are fitted by the sum of signal and background 
functions, described by a Gaussian and a product of an exponential and a polynomial, 
respectively. In case of the D* sample a second Gaussian is used to describe the reflection 
of_D° — > Kirn decay, where the ir° meson is not observed. In the D° sample this 
reflection is not visible. The total number of D° mesons is about 3,800 and 13,800 in the 
D* and the D° samples, respectively. 

Finally, the determination of (Ag/g) x requires knowledge of the analysing power 
Oll- Since only one D° is measured, the partonic kinematics cannot be fully reconstructed 
and Oll cannot be calculated on an event-by-event basis. A kinematic factor which is 
approximately equal to the depolarisation factor D is factored out and the remaining 
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Table 1: Systematic error contributions to (Ag/g) x . 



source 




source 




False asymmetry 


0.09 


Beam polarization Pb 


0.02 


Fitting 


0.09 


Target polarization P t 


0.02 


Binning 


0.04 


Dilution factor / 


0.02 


MC parameters 


0.05 







Total error 0.15 



part of oll is parameterized in terms of measured kinematic variables. This is done using 
a neural network trained on a Monte Carlo (MC) sample for D* mesons. The correlation 
between the generated and the reconstructed oll is 82% (see Fig. [2]). The sample was 
generated with AROMA [19] in leading order QCD and the events were processed by 
GEANT to simulate the response of the detector and finally reconstructed like real events. 
The scale, /z 2 , used in the MC simulation was determined by the mass of the produced 
charm quark pair and is sufficiently large to justify the perturbative approach. 

The gluon polarisation (Ag/g) x is determined for each of the 29 weeks of data taking 
with a standard least square minimisation procedure taking into account the statistical 
correlation between events weighted by ws and by As in our inclusive analysis [I] the 
event weighting reduces the statistical error. The final value for (Ag/g) x is the weighted 
mean of the above results. The resulting As is consistent with zero. Note that we measure 
Ag/g in a given range of x. Provided that Ag/g(x) is weakly dependent on x in the range 
covered, this method gives a measurement of Ag/g((x)), where (x) is calculated with 
signal weights. The above assumption is supported by the results of our QCD analysis [I]. 

The major contributions to the systematic error are listed in Table [TJ The contribu- 
tions from Pb ; P and / are as discussed in [JJ. To study the influence of false asymmetries 
the final samples from Fig. 1 were subdivided into two samples using criteria related to the 
experimental apparatus, e.g. kaons going to the upper or to the lower spectrometer parts. 
The resulting asymmetries were found to be compatible within their statistical accuracy, 
thus no false asymmetries were observed. The upper limit of the contribution to the sys- 
tematic error was also estimated from the dispersion of the values for (Ag/g) x and As 
for the various data taking week. Assuming that possible detector instabilities are similar 
for background and signal events and applying a method as in [4] leads to a conservative 
limit of 0.09 for both decay channels. The fit for the signal purity determination was per- 
formed with different background parameterizations, different binnings for the invariant 
mass spectra and varying the constraints for some of the fit parameters. The resulting 
dispersion of the (Ag/g) x values was 0.09 for both channels. The (Ag/g) x calculations 
were repeated with several sets of binning in /P&gll an d the dispersion of the results was 
0.04 for both channels. Other contributions like radiative corrections and event migration 
between target cells are negligible. All these studies were done independently for the D° 
and D* samples and result in very similar values for all the contributions. To estimate the 
influence of the simulation parameters, i.e. charm quark mass (1.3 GeV/c 2 to 1.6 GeV/c 2 ) 
and parton distribution functions, MC samples with different parameter sets were gener- 
ated and oll was recalculated. The dispersion of the resulting (Ag/g) x for the combined 
D° and D* sample was 0.05. In addition, it was checked that the parameterization of oll is 
valid for the D° and the D* sample. The resolved photon contribution to the open charm 
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Figure 3: Compilation of the {Ag/g) x measurements from open charm and high hadron pair 
production by COMPASS [13], SMC [H] and HERMES pi] as a function of x. Horizontal bars 
mark the range in x for each measurement, vertical ones give the statistical precision and the 
total errors (if available). The open charm measurement is at a scale of about 13 (GeV/c) 2 , other 
measurements at 3 (GeV/c) 2 . The curves display parameterizations from a NLO QCD analysis 
in the MS scheme at 3 (GeV/c) 2 , [4J: fits with AG > (broken line) and with AG < (dotted 
line). 



production via gluon-gluon fusion has been estimated with the RAPGAP generator [20] 
and found to be negligible in our kinematic range. 

The values obtained for (Ag/g) x are 0.53 ± 0.75(stat) for the D° and —1.01 ± 
0.55(stat) for the D* sample. These two values agree within 1.7 standard deviations and 
their weighted mean is 




= -0.47 ± 0.44(stat) ± 0.15(syst) 



at a value of (x) = 0.1 l-o.ol an d a scale (/i 2 ) ~ 13 (GeV/c) 2 . All the contributions to the 
systematic error in Table Q] were added in quadrature and conservatively assumed to be 
fully correlated in the two samples. 

In Fig. [3] the above result is compared to other measurements of (Ag/g) x and to 
parameterizations from a QCD analysis of the structure function data 0]. It is in good 
agreement with previous measurements favouring small values of (Ag/g) . Note that 
the scale here is /i 2 ~ 13 (GeV/c) 2 while all other points and the curves are given at 
fi 2 « 3 (GeV/c) 2 . 
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In summary, we have performed the first determination of (Ag/g) x from a measure- 
ment of the cross section asymmetry for D° meson production. In the analysis photon- 
gluon fusion in LO QCD was assumed to be the underlying production mechanism for 
open charm production. The resulting value of (Ag/g) x is compatible with our previous 
result from high hadron pairs but is less model dependent. 

We acknowledge the support of the CERN management and staff, the special effort 
of CEA/Saclay for the target magnet project, as well as the skills and efforts of the 
technicians of the collaborating institutes. 
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